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Abstract. We describe some of the recent progress in the calculation of ther- 
modynamic quantities in QCD at high temperatures and densities by weak- 
coupling techniques and extrapolation to realistic coupling strength. We argue 
that a (mostly) weakly coupled quark-gluon plasma at temperatures only a few 
times the transition temperature is not yet ruled out by the observed fast appar- 
ent thermalization at RHIC, as nonabelian plasma instabilities might provide 
an efficient mechanism for fast isotropization even in a collisionless plasma. 
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1. Introduction 

To approach the thermodynamics of the strong interactions by weak-coupling tech- 
niques [EElEl might appear to be a hopeless enterprise, at best of academic interest 
or applicable only to such high temperatures and densities that it is completely irrel- 
evant for the hunt of the quark-gluon plasma in heavy-ion physics. However, while 
the physics of the deconfmement transition certainly cannot be captured by Feyn- 
man diagrams of quarks and gluons, it may actually be the case that extrapolations 
from asymptotically high temperatures already begin to make sense quantitatively 
at temperatures only a few times the dcconfincmcnt temperature, and I shall present 
some evidence for this in what follows. It is true that strict perturbation theory fails 
spectacularly, even at temperatures as high as 10 5 times the deconfmement tem- 
perature, where perturbative QCD should work without difficulty. But this failure 
has little to do with nonperturbative properties of nonabelian gauge theories — it 
even occurs in such trivial theories as scalar field theory in the large- N limit, where 
thermal quasiparticles are free of interactions. Indeed, this failure can be repaired 
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to some extent (at least for bulk thermodynamics), as I shall describe below, and 
the results from comparison with lattice results are quite encouraging. 

It may be argued that the partial successes of weak-coupling techniques at 
temperatures a few times the deconfinement temperature do not prove anything, and 
that the quark-gluon plasma at such temperatures is essentially strongly coupled. 
This is by now certainly the opinion of the majority since the strongly coupled 
quark-gluon plasma has been heralded as one of the major discoveries at RHIC. 
The evidence quoted for this is the fast apparent thermalization and quick onset 
of hydrodynamic behaviour, which seems incompatible with existing perturbative 
predictions [0J |SJ 03 El El IH] ■ However, it is not excluded that weak-coupling 
physics could explain the fast apparent thermalization. Arnold et al. [El El EH 
have recently pointed out that there is only evidence for fast isotropization, and 
the latter may occur also in a weakly coupled plasma through nonabelian plasma 
instabilities as advocated since long by Mrowczyriski [ 1151 1161 I17| . Most recently, 
Romatschke, Strickland, and myself [ I18| have produced supportive evidence for this 
picture by numerical simulations of nonabelian plasma instabilities in the nonlinear 
hard- loop approximation [ I19j . The latter is a weak-coupling framework, but it 
nevertheless produces nonperturbative effects in the form of nonperturbatively large 
fields which can efficiently wipe out anisotropies in a collisionless plasma [E]- Fig-D 
shows the exponential growth of 1+1-dimensional nonabelian plasma instabilities 
(constant modes in the transverse plane) which continues into the nonlinear regime 
so that it can only be stopped when they begin to affect the trajectories of hard 
particles. 




(a) (b) 

Fig. 1. (a) The exponential growth of nonabelian plasma instabilities as measured 
by the rms value of the induced current density j (full line) and a measure of 
local abelianization (C) in the numerical simulation of Ref. [ I18|: (b) the total 
energy transferred from hard modes into the soft unstable modes (full line) and the 
energies in the various transverse (longitudinal) chromoelectric (-magnetic) field 
components. 
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2. Resumming perturbative results at high temperature 

However, even if this picture is not complete, and if there are important nonper- 
turbative effects, it is certainly worth knowing how much of the physics of the 
quark-gluon plasma can already be described by weak-coupling techniques. The 
answer will certainly depend on the quantity under consideration, and the easiest 
quantities should in fact be those related to the equation of state, as the thermody- 
namic potential should be dominated by contributions from hard modes, for which 
the effective coupling should be smallest. Indeed, comparing perturbative results 
with lattice results, the leading interaction contributions to the thermodynamic 
potential oc g 2 are doing quite well in reproducing the measured deviation from 
ideal-gas behaviour for T > T c , see FigEl 
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Fig. 2. Strictly perturbative results for the thermal pressure of pure-glue QCD as 
a function of T/T c (assuming T c /A^g = 1.14). The various gray bands bounded by 
differently dashed lines show the perturbative results to order g 2 , <? 3 , g 4 , and g 5 , 
using a 2-loop running coupling with MS renormalization point /2ms varied between 
7rT and AirT . The thick dark-grey line shows the continuum-extrapolated lattice 
results from reference [EI]; the lighter one behind that of a lattice calculation using 
an RG- improved action [I22|. 




Problems arise from the next terms in the perturbative expansion, where collec- 
tive phenomena such as Debye screening come into play [^. Up to and including 
order g 5 the thermal pressure of a quark-gluon plasma with Nf massless quarks 
reads reads [E3 EH EH] 



P = ^T 4 { (1 + %N f ) - f (1 + &N f ) £ + 30 [(1 + \N f ) 
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where a s = g 2 /(47r) and /2ms is the renormalization scale parameter of the MS 
scheme. The higher-order terms involve softer scales, but the problem is not that 
the effective coupling of these modes is larger. Even when the running coupling 
is chosen as that appropriate for hard modes, the next term, which is of order g 3 , 
comes with a coefficient such that only at temperatures much larger than 10 5 T C 
there is apparent convergence of the perturbative series. At the same time, the 
dependence on the renormalization scale increases instead of decreasing, and this 
does not get better when higher terms in the perturbative series are added in, quite 
to the contrary. 

However, as I have already mentioned, this complete failure of the perturbative 
expansion is not related to specifically nonabelian physics. Rather, it can be found 
in almost any thermal perturbation theory, such as in the rather trivial large- N limit 
of scalar field theories. The latter can be "solved" exactly, yielding a rather boring 
thermodynamic potential as a function of the coupling, and yet the perturbative 
approximations oscillate wildly and seem incapable of describing the full result 
except at tiny coupling [ I26| . In this example clearly no new physics needs to 
be incorporated, but thermal perturbation theories should be reorganized. One 
proposal for doing so is called "screened perturbation theory" [|23 EE], and has 
been generalized to QCD by Andersen, Braaten, and Strickland [ 1291 I30| . An 
alternative proposal, by Blaizot, Iancu, and myself, is based on an expression for 
the entropy density that can be obtained from a ^-derivable two-loop approximation 
[ |3*T1 152] (see also These approaches indeed succeed in taming the plasmon 

term ~ g 3 that spoils the apparent convergence of strict perturbative expansions in 
g. Moreover, in Refs. [ 1311 132} the lattice results for the thermodynamic potential for 
T > 3T C were quite well reproduced by transforming the leading-order interactions 
into so-called hard-thermal- loop [ 134) quasiparticle properties (which include Debye 
screening). 

The potential for improvements of the perturbative results through partial re- 
summations can also be seen in the a priori strictly perturbative framework of 
dimensional reduction. Using these techniques, the perturbative expansion of the 
QCD thermodynamic potential has by now been carried out up to and including 
order g 6 ln(g). The result is [ I35| 



Including this term in a strictly perturbative expansion only worsens the already 
disastrous apparent convergence. However, the results Q and J2J each involve 
contributions from quite different sources. One contribution is from hard modes, 
and this behaves decently as concerns apparent convergence and renormalization 
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scheme dependence. The problematic part is from the soft sector, which is governed 
by an effective three-dimensional Yang-Mills theory with the chromoelectric degrees 
of freedom turned into adjoint scalars with mass wle- Its contribution to four-loop 
order is 
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where c is a constant that is inherently nonperturbative and needs 3-d lattice cal- 
culation to be determined. 
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Fig. 3. (a) Three-loop pressure in pure-glue QCD with unexpanded effective-field- 
theory parameters when /2ms is varied between ttT and 4rrT (medium- gray band). 
The broad light-gray band underneath is the strictly perturbative result to order g 5 
with the same scale variations. The full line gives the result upon extremalization 
(PMS) with respect to p,MS (which does not have solutions below ~ 1.3T C ); the 
dash-dotted line corresponds to fastest apparent convergence (FAC) in m%. (b) The 
same extended to four-loop order by including the recently determined g & \n.(l/g) 
contribution of together with three values for the undetermined constant 5 in 
[g 6 ln(l/<7) + 5]. The broad light-gray band underneath is the strictly perturbative 
result to order g 6 corresponding to the central value S = 1/3, which has a larger 
scale dependence than the order g 5 result; the untruncated results on the other hand 
show rather small scale dependence. The full line gives the untruncated result with 
5=1/3 and /2ms fixed by PMS (which does not have solutions below ~ 1.9T C ); the 
dash-dotted line corresponds to fastest apparent convergence (FAC) in m E . (Taken 
from [Eg) 



The effective field theory parameters wie and gE can be calculated by pertur- 
bative matching, and if their dependence on g is expanded out and truncated at a 
given perturbative order, the bad convergence properties discussed above arise. 

However, the simple prescription of keeping effective field theory parameters 
unexpanded in g leads to considerable improvements [I3fi|. which are most striking 
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Fig. 4. (a,) Comparison of the exact result for the pressure of large- Nf QCD 
with perturbation theory to order g 2 , g 3 , g 4 , and g 5 , respectively, at different 
renormalization points (taken from [Ej); (b) same exact result, but as a func- 
tion of 5cff(^MS = tT) and, alternatively, log 10 ( Al/7tT) . The purple dashed line 
is the perturbative result when the latter is evaluated with renormalization scale 
A*ms = Mfac = 7re 1 / 2_7 T; the blue dash-dotted lines include the numerically deter- 
mined coefficient to order g® s (with its estimated error) at the same renormalization 
scale. The result marked u g^ s = g®g v corresponds to choosing /2ms such that the 
order-<^ ff coefficient vanishes and retaining all higher-order terms contained in the 
plasmon term ex m 3 E . The (tiny) red band appearing around the exact result at 
large coupling displays the effect of varying the cut-off from 50% to 70% of the 
Landau-pole scale Al (taken from [EHj)- 

when the term of order g 6 ln(g) is included and the undetermined constant under 
this log in (J2J) is fixed. This can be seen in Fig. El which even suggests a definite pre- 
diction for this constant under the log, which hopefully will one day be determined 
by a combination of analytical and lattice techniques. 

A test bed for resummations of perturbative results which has many of the spe- 
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cial features of gauge theories at finite temperature and/or chemical potential is the 
large flavour-number limit of QED or QCD, where the thermal pressure can be calcu- 
lated exactly [^7\ 12H] • Fig-EK shows the lack of convergence of a strict expansion 
of the thermodynamic pressure in powers of g and the large renormalization-point 
dependences. Fig. on the other hand shows that the above-mentioned prescrip- 
tion of keeping unexpanded the effective field theory parameters in the dimensional 
reduction result to order g 6 reproduces the exact result surprisingly well for all 
couplings for which the large- Nf theory is still well-defined (at the highest coupling 
shown the Landau pole of this theory is being approached). 

3. Finite quark chemical potential 

In recent years there has also been substantial progress in exploring the equation 
of state of QCD at high temperature and nonvanishing quark chemical potential 
on the lattice [0U[ ^fl 021- These results do not yet have the quality of results at 
zero chemical potential, and no reliable continuum extrapolations have been carried 
out yet. But it is plausible that these results are already approximating continuum 
results with errors of less than about 20%. Comparing HTL-quasiparticle resum- 
mation results [031 \SM or improved results from dimensional reduction techniques 
extended to finite chemical potential EE3 E±] with these lattice data is in fact 
most encouraging. Even at the rather low temperatures that are involved, the ana- 
lytical weak-coupling results reproduce the equation of state within expected errors, 
see Fig. □ [EH]. 

There has been for some time a discrepancy between lattice results on off- 
diagonal quark-number susceptibilities (i.e. derivatives with respect to chemical 
potentials for different quark flavours) and the perturbative result, whose leading- 
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Fig. 5. The difference AP = P(T, /j,)-P(T, 0) divided by T 4 using the unexpanded 
three-loop result from dimensional reduction of Ref. [E| for n/T = 0.2, . . . , 1.0 
(bottom to top). Shaded areas correspond to a variation of /2ms around the FAC-m 
choice by a factor of 2; dashed and dash-dotted lines correspond to the FAC-g and 
FAC-m results, respectively. Also included are the recent lattice data of Ref. [ 141) 
(not yet continuum-extrapolated!) assuming T^ =0 — 0.49 Aqcd- 
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Fig. 6. Complete entropy density in the large- iV/ limit for the three values 
3cff(MMS = 2/i) = 1,2,3 (heavy dots), compared with the full HDL result (solid 
line) and the low-temperature series up to and including the T 3 In T contributions 

[Eg. 

order term turns out to be of order g 6 \n(g) [02]. The authors of Refs. [EI3 1ST) 
have obtained results in the deconfincd phase that were far below those predicted by 
perturbation theory, and have interpreted this as new evidence for nonperturbative 
physics and the failure of weak-coupling methods. Most recent lattice results [ 
1521 1531 151| have disproved the previous ones, and there is now agreement with the 
perturbative estimate at T > 2T C . 

With lattice methods, it is possible to explore only moderate quark chemical 
potentials at high temperature. Using weak-coupling techniques and extrapolations 
down from asymptotically large chemical potential, much progress has been made 
in recent years in unravelling the rich physics of QCD at low temperatures and 
high quark chemical potential, which seems to involve numerous different phases of 
colour superconductivity [ EHI EB E3 • 

Above the critical temperature of colour superconductivity, which is expected 
to be somewhere between 6 and 60 MeV [ EH , the thermodynamic behaviour is 
in fact still fundamentally different from that at high temperatures. The only 
weakly screened chromomagnetic interactions lead to non-Fermi-liquid behaviour [ 
EH1 EI3 EH1 E2 En] which among other effects leads to an anomalous specific heat 
when T becomes comparable to or smaller than gfi. Whereas this effect has been 
discussed in QED a long time ago [ I64 | with the conclusion that this effect is probably 
unobservably small, in QCD the corresponding effect is orders of magnitudes larger, 
since the strong structure constant times the number of gluons is much larger than 
the fine structure constant times one. Until recently, however, the corresponding 
resummed perturbation series has been unavailable except for the coefficient in 
front of a leading logarithmic term [ 65 . This has been completed in [ ESI EH j 
leading to a low-temperature entropy with a perturbative expansion involving logs 
and fractional powers of the temperature of the form 
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+ 1.041 3 f^J 3 + 0(T 3 InT), T « 5M . (4) 

Using hard-dense- loop (HDL) resummations [E3 EH EH], quantitative results 
are also available for T ~ [[HS] and these can be tested against the solvable case 
of l&rge-Nf QCD, as shown in Fig. Physically interesting applications of these 
results are presumably outside the reach of heavy-ion colliders. But they may be 
of interest to neutron star physics, namely the cooling behaviour of young neutron 
stars if they have a normal-conducting quark-matter component. Both the specific 
heat and the neutrino emissitivity contain large logarithms from non-Fermi-liquid 
behaviour, with important deviations [ 1701 171j from naive lowest-order perturbation 
theory [172], 

4. Conclusions 

To summarize, there has been quite some recent progress in evaluating by weak- 
coupling techniques the thermodynamics of a deconfined quark-gluon plasma. At 
zero chemical potential, the calculation has been pushed to order g 6 log(g), and at 
large chemical potential analytical results exist for the non-Fermi-liquid behaviour 
of ungapped quark matter and also of colour superconductors. At zero to small 
chemical potential, comparisons to lattice results are possible and for T > 2T C seem 
to validate the outcome of resummations that are necessary to overcome the poor 
apparent convergence of strict perturbation theory. 
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